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Abstract 

We show that the exchange statistics have consequences in ID systems with 
compact topology, contrary to the common opinion that exchange statistics is ar- 
bitrary in ID. As examples of non-trivial statistical behavior we exactly calculate 
the partition function and correlators for systems of free q-particles on compacti- 
fied chains using functional integral techniques and the supersymmetric trick. In 
particular we consider a spin 1/2 XY-chain with periodic boundary conditions that 
corresponds to the case of q=-l. 



Last decade there has been considerable interest in various deformations of quantum 
statistics which interpolate between the Bose- and Fermi- cases. Such deformations can be 
separeted into the deformations of exchange phase which appears under the permutation 
of particles (exchange statistics with braiding [|I| or without it 0) and the deformations 
of Pauli principle (so called exclusion statistics 0). While the last can be defined in 
any space dimension the possibility to introduce non-trivial exchange statistics crucially 
depends on the dimension. 

Moreover there is a common opinion that the exchange statistics is something arbitrary 
in ID. It is due to the fact that in ID for the particles with hard-core condition on a 
finite interval there no real permutation is allowed and, hence, the exchange phase never 
plays a role. The situation changes if the interval is reduced to a circle with periodic 
boundary conditions. This leads to a non-trivial loop which gives new possibilities to 
permute particles via the 'glued' boundary. Then the effects of the exchange statistics 
immediately exhibit themselves as we will demonstrate in this Letter. However we should 
add that we deal with lattice systems because the effects of the loops are negligible in 
continuous limite. 

The present calculations use the functional integral formalism. This is not surprising 
if we remember that the formal form of the functional integral representation for the 
thermodynamical quantities are essentially the same for both Fermi- and Bose- statistics 
and the difference is reflected only in the nature of fields and boundary conditions We 
will now show that the statistics can be taking into account by changing the boundary 
conditions both in imaginary time and space which lead to the transmutation of distri- 
bution functions. In the framework of our model we calculate exactly such distribution 
functions and, moreover, the result is simple. This hints at the form of multidimensional 
analogues (see discussion after Eq.(|)). 

The point is illustrated by the example of free q-fermions, i.e. the particles with 
deformed exchange statistics defined by the algebra of commutation relations for the 
creation (annihilation) operators af, (oj) at the i-th site of the compact chain with length 
M. The relations are split into two parts: the first describes the fermionic algebra on 
each site: 

afai + ttittf = 1 , = , (ai)+ = af , (1) 

and the second gives the exchange statistics with deformation parameter 
r > 1: 

ttiaj = —qajai , afaj = —q^^ajaf , (2) 

where the ordering 1 < i < j < M is assumed. 

Presently we will consider a gas of free q-fermions with the Hamiltonian: 

M 

H = ^{Ai^i-iaf Qi-i + l/2Biafai) + qAiMafau + l/2Bialai + H.c. (3) 

1=2 

One of the main purposes of the Letter is to examine the consequences of g-deformed 
exchange statistics for the thermodynamical quantities on the basis of the Hamiltonian 
(^. It turns out that the partition function and correlators for the model can be calculated 
exactly. So we are going to do for the exchange statistics the same that was done by Ha 
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for the exclusion statistics: to examine statistics from the point of view of exact solutions. 
It should shed light on the whole problem. 

First of all let us describe and discuss the main results for the thermodynamical quan- 
tities associated with the Hamiltonian (^) and then give the outline of the calculations. 
The starting point for this is the exact expression for the partition function of the sys- 
tem, Z = Tt exp{—f3{H — fiN)), at temperature 1/(3 and chemical potential /i, where 
I aftti is the number operator for q-particles: 

-| n—l n—1 

z = -J2J2 1''' 11(1 + gPe'^(/^-'('=r))) (4) 

^ 1=0 p=0 fc™ 

Here the notation e'(/c[") was introduced for the kf^-th eigenvalue of the hopping matrix 
\\Aij\ \ on the circle with g'-periodic boundary conditions (which corresponds to the sub- 
stitution of q^'AiM instead of Aim in the hopping matrix). For the homogeneous chain 
[Aij = A) these eigenvalues have the form: 

e\kr) =B + 2Acos{kr) , = ^(m - Ir/n) (5) 

with m = 0, . . . , M — 1. Using expression (H) we immediately obtain the representation 
for the distribution function /i) of the q-particles in terms of transmuted Fermi-Dirac 
(or Bose-Einstein) functions fg{x) = 1/(1 + g~^exp(a;)): 

So we see that in the Eqn.(|^) q-deformed Fermi-Dirac functions, with all allowed powers 
of the deformation parameter q, make contributions. Such functions naturally interpolate 
between fermionic distribution function {q = 1) and bosonic distribution function (g = 
— 1). Moreover they pick up all g-periodic boundary conditions which also reflects q- 
exchange factor gained due to the permutation of particles via the loop. 

The note about all g-periodic boundary conditions let us to look at the problem from 
the different position. Indeed, from the original form of the Hamiltonian (^, which is 
quadratic on the creation (annihilation) operators, it is not difficult to see that, using 
the appropriate diagrammic technique (for example in the g-field technique 0), all con- 
tributions to the thermodynamical quantities will be given by loop diagrams with some 
exchange factors. And this does not depend on the space dimension. On the other hand, 
the examination of the expression (^) for the ID case may prompt an ansatz for the 
partition function in nD after the summation of such contributions for different loops: 
the sums over all loops in the configuration space, with all allowed g-periodic boundary 
conditions for them and all possible g-deformations of Fermi-Dirac constructions. 

Let us turn our attention to the particular case g = — 1. In this case the Hamiltonian 
(^ can be realized as a Hamiltonian of the spin 1/2 compact XX-chain in a magnetic 
field B (for the sake of simplicity we deal with the homogeneous case): 

H = Efi2 {A,s^s^ + Aysfs]-Bs^,)-B,sl 

+A,sfsl, + Ayslsl + MB/2 ^ ' 
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if we identify = Ay and the creation (annihilation) operators with spin upper (lower) 
operators are usual: 



a.: 



sf + isl 



After the introduction of anisotropy Ax{y) = 2(A ± 7) Hamiltonian (|^ becomes the 
Hamiltonian of a spin 1/2 compact XY-chain in the magnetic field B. As we will show 
below, we may obtain the exact formulae for the thermodynamical quantities for the 
Hamiltonian of XY-model in compact lattice topology (in contrast to previous works 
where such quantities were calculated in thermodynamical limite |10|, |ll[). Now 

we only state them. The partition function of the model Z{f3, B) contains four terms, two 
of which have fermionic nature and the others have bosonic ones: 



z = -(z; + z7 + i/z+-i/z, 

Here Z^[Z^) are fermionic (bosonic) partition functions 



(9) 



7± 



-PBM/2 



^± ^ g/3BM/2 -Qg-/3£;(fc±)/2 

for systems with the energy spectra 

E{k^) = (fi + 2Acos(A;^)). 



/(I 



(27sin(fc±))2 
(5 + 2Acos(A;±))2 



and antiperiodic (periodic) boundary conditions: 



k- 



2tt 



m + 1/2) , k' 



2tt 
M 



m m = 0,...,M — 1. 



We also calculated correlators for various spin components. Since the results are very 
cumbersome, they will be reported elsewhere. 

Now let us outline the calculational procedure which leads to expressions (0),(|) and 
then give the modification for the case of XY-model (0). To this end we apply to the cre- 
ation (annihilation) operators af, (a^) some analog of the Jordan- Wigner transformation 
just as it was performed in 1961 by E.Lieb, D.Mattis and T.Schults in Ref. |0] for the case 
of the Paulion chain (q=-l): 



cfq Sfc=i4cfe 



(10) 



Where the above operators cf, (q) are respectively creation (annihilation) operators of 
spinless fermions at the i-th site. The operator (|^) is cast in the form: 



M 



H = ^(A,i-ic+Q_i + l/2Bictci) + AiMq^cfcM + l/2Bi4ci + H.c. 



'IV. 



1=2 
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where N = X^i^fi cfci is the number operator of the fermions. To proceed we introduce 
the set of operators {-P;}^o^ which are projection operators on the subspaces of /(modn) 
particles: 



n~l 



- E Q'^^'Q^'' ■ (12) 



Making use of these projection operators with the obvious property J2?=o Pi = I ^ve then 
obtain the following representation for the partition function Z = Ti exp{—j3{H — nN)): 

n-l 

Z = E Tr expi-PiH - ^N))Pi (13) 

1=0 

that allows to replace in each term of the sum (|T3|) the multiplier in the Hamiltonian 
( pUD by the numerical constant qK Moreover due to the commutation of the operators H 
and N each p-th term of the expression (|T^) for the projection operators can be taken 
into account by a shift of the chemical potential /i to the value fj.^^ = A* + 27ripr/np. This 
leads to the representation for the partition function: 



-1 n—l n—1 

Z = - E E ^''"Tr exp(-/3(if(') - fi^^^N)) (14) 
;=0 p=0 

where if*^'-* are g'-periodical square fermionic Hamiltonians. After the making of use of 
the standard results for the partition function of a free fermion system we regain to the 
Eq.(^). In the same way various averages and correlation functions can be obtained. For 
example: 

-tMB n-l n-l 

(ai(r)a+^,(0)) = E E ^^'^D^t Cj(r)(Cj(r))^i^, 

1=0 p=0 

where the matrix C is defined by the relation 

+£ Efcr e(--^)(^'('=r)-/^)+^'=rK-™")(i + (g-1 _ i)0(L - m')) ^ ' 

(here we used lattice 9 symbol which is continuous from right). It is not difficult to check 
that this average reproduces the formula (^ for the distribution function if r, L — > 0. 

By the same way we calculate the generating functional (e"'^^^-') with Q{L^ — Efc=i c^^Cfc. 
Using the expression (^) for the partition function Z we have 

X*^ / — nZ ^1=0 ^p=0 y 

Up to this point we did not use the functional integral approach. However it allows 
us to look at the problem from another viewpoint and, furthermore, is closely connected 
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with the modification of the calculations for the case of XY-chain. Indeed, each term of 



the sum (|T4|) can be represented in the form of a function integral [Q: 
Tr exp(-/3(i/(') - /x^^'^A^)) = J D^{T)D^{T)e^'r> 

over the Grassmann fields ^{t),^{t) with antiperiodic boundary conditions ^(/?) = — ^(0), 
= —^(0) and the action 

Sip= !^ dri^^^ - H^%U) + m + !^ dr'-^U ^^^^ 
= Si + 6Sp 

— — ■ 27Tpr ■ 2-Kpr 

Under the change of variables ^(r) ^(r)e~*^^^, ^(r) — >■ ^(r)e*^^^ the last term of the 
previous expression, 5Sp, disappears, that, however, is compensated by the changing of 
the boundary conditions for the fields of integration: 

m = -q-'m , m = -i^m , (is) 

(a procedure of such type was used in Ref. [|T^ to construct the Feynman diagram 
technique for spin systems). As a result we reach the following functional integral repre- 
sentation of the partition function: 

n— 1 n— 1 „ 

^ = E E q''' / Dar)Dar)e'^/n (19) 
p=0 1=0 •' 

provided with the g^-periodic boundary conditions ([18|). In this form the equation turns 
out to be valid for the case of XY-model for spin 1/2, although the justification needs to 
attract other ideas. 

Let's clarify the problem arising in the case of XY-model. Indeed, after the Jordan- 
Wigner transformation (|10D the Hamiltonian (|^) becomes the operator: 

H = E-4(Ac+Q_i + l/25,c+Q + 7c+cti) + l/2B,ctc, 

~Ai-lf ctcM - li-lf ctct, + H.c. ^^^^ 

which differs from the Hamiltonian (0) by additional terms with the anisotropy param- 
eter 7. They cause the Hamiltonian to no longer commute with the number operator 
for particles A^, that does not allow us to use device of shifting the chemical potential. 
However the Hamiltonian still commutes with the operator r = (—1)^. This leads to the 
formula, which is analogous to the Eqn.(0): 

Z = |( Tre-/5(^^<+'-MJV) + Tre-^(^'"'-'^^) + 

rYj.^^-l3iH(+)-f,N) _ rp^^g-/3(/f(-)-/.iV)^ 

where if*^^^ coincide with the Hamiltonian (|20|) after the substitution r = (—1)^ ±1. 
Although the two last terms cannot be calculated by changing of chemical potential, 
taking into account r in them is also achieved by changing of boundary conditions ([T8|). 
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This was demonstrated in the Ref. |13| into the framework of supersymmetric quantum 
mechanics, where the operator r and Tr re^^^ play the role of supersymmetric involution 
and supersymmetric Witten index of the Hamiltonian H. That is why we term this the 

supersymmetric trick. 

All this returns us to Eg . (|T9D with g = — 1, which after the use of the Bogoliubov 
transformation for the variables of integration leads to the formula (P). We want to add 
that the same trick was used to calculate correlation functions. 

The Letter concentrates essentially on the principal questions concerning exchange 
statistics. However the formulae obtained in the Letter can have various physical appli- 
cations. Let's briefly mention them. They describe the compact finite chain in contrast 
to the classical papers [0, § ^ and more recent ones |T^, ^ where thermodynamical 
quantities were calculated in the thermodynamical limit. In this limit the statistical ef- 
fects, which are proportional to the inverse size of the system, 1/M, are negligible and the 
formulae of the paper are transformed into the known results. However recently a lot of 
attention was attracted by the theory of so-called J-aggregates, i.e. molecular aggregates 



with an unusually sharp absorption band ([|1J],[|T5[ and Refs. wherein) where the advances 
were connected with the use of exact results on ID chains. Frenkel excitons in such long 
molecules obey the Paulionic commutation relations [1^, or, in a more general case (if we 
will take into account retardation effects ^^), q-particle commutation relations. We hope 
that theformulae obtained will be used in the theory of the nonlinear response of long 
cyclic molecules where a thermodynamical limit is not appropirate. On the other hand 
XY-model recently arose in models of adsorption processes with diffusional relaxation 



1^. Here the forlulae of this Letter can also find applications to describe the compact 
case. The same one can be said about the application of results to the theory of defectons 



where the statictics of defectons is exactly Paulionic statistics! 19 



In conclusion, we exactly calculated the partition functions and correlators for systems 
of free particles with q-deformed exchange statistics and the spin 1/2 XY-model on a 
compact chain with periodic boundary conditions. We demonstrated that the deformation 
of the statistics can be taken into account by changing the boundary conditions for the 
fields of integration in the functional integral framework that lead to the transmutation of 
distribution functions. In such a way deformed Fermi-Dirac functions appeared. We also 
stated a conjecture about the form of thermodynamical quantities in higher dimensions. 
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Grants N 94-02-03712 and N 95-01-00548 and partially (K.I.) by the UK EPSRC under 
Grant number GR/J35221, Euler stipend of German Mathematical Society and grant 
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